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Abstract 

We study the decays of the lightest top squark in supergravity models with and 
without R-parity. Using the simplest model with an effective explicit bilinear break- 
ing of R-parity and radiative electroweak symmetry breaking we show that, below 
the threshold for decays into charginos t\ — ► b\\, the lightest stop decays mainly 
into third generation fermions, t\ — > br instead of the R-parity conserving mode 
ti ~~ > c X?j even for tiny tau-neutrino mass values. Moreover we show that, even 
above the threshold for decays into charginos, the decay t± — > br may be dominant. 
We study the role played by the universality of the boundary conditions on the 
soft supersymmetry breaking terms. This new decay mode t± — ► br as well as the 
cascades originated by the conventional t\ — > cxi decay followed by the R-parity 
violating neutralino decays can provide new signatures for stop production at LEP 
and the Tevatron. 



PACS: 14.80.Ly, 11.30.Pb, 11.30.Fs, 12.10.Dm, 12.10.Kt, 12.60.Jv 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) |jj or its minimal supergravity 
(SUGRA) version |2] are by far the most well studied realizations of supersymmetry. 
However, neither gauge invariance nor supersymmetry requires the conservation of R- 
parity. Indeed, there is considerable theoretical and phenomenological interest in studying 
possible implications of alternative scenarios @ in which R-parity is broken [|], |5], || |7j. 
This is especially so considering the fact that it provides an appealing joint explanation 
of the solar and atmospheric neutrino anomalies which has, in addition, the virtue of 
being testable at future accelerators like the LHC The violation of R-parity could 
arise explicitly as a residual effect of some larger unified theory |J, or spontaneously, 
through nonzero vacuum expectation values (vev's) for scalar neutrinos ||, || 0. In 
realistic spontaneous R-parity breaking models there is an SU(2) <g>U(l) singlet sneutrino 



vev characterizing the scale of R-parity violation |10|, [II], [1^, [13[ which is set by the 
supersymmetry breaking scale. 

There are two generic cases of spontaneous R-parity breaking models to consider. In 
the absence of any additional gauge symmetry, these models lead to the existence of a 
physical massless Nambu-Goldstone boson, called majoron (J) which is the lightest SUSY 
particle, massless and therefore stable. It plays an important role in making these models 
fully consistent with astrophysics and cosmology || if one wishes to contemplate the case 
of large breaking scales and heavy tau neutrino. If the majoron acquires a small mass due 
to explicit breaking effects at the Planck scale then it may decay into electron and muon 
neutrinos or photons, on very large time scales of cosmo logical interest, playing a possible 
role as unstable dark matter |TJ]]. Alternatively, if lepton number is part of the gauge 
symmetry and R-parity is spontaneously broken then there is an additional gauge boson 
which gets mass via the Higgs mechanism, and there is no physical Goldstone boson |Tj| . 
As in the standard case in R-parity breaking models the lightest SUSY particle (LSP) 
is in general a neutralino. However, it now decays mostly into visible states, therefore 
diluting the missing momentum signal and bringing in increased multiplicity events which 
arise mainly from three-body decays such as 

X°x - ffr, (1) 
where / denotes a charged fermion. The neutralino also has the invisible decay mode 

X? - 3i/. (2) 

as well as 

xl^vJ, (3) 



in the case the breaking of R-parity is spontaneous [10], pTT| . This last decay conserves 
R-parity since the majoron has a large R-odd singlet sneutrino component. 

If R-parity is broken then supersymmetric (SUSY) particles need not be produced 
only in pairs, and the lightest of them could decay. The effects of R-parity violation can 
be large enough to be experimentally observable. 

In this paper we focus on the decay modes of the lightest top squark in supergrav- 
ity models where supersymmetry is realized with R-parity violation. In such models the 
lightest stop could even be the lightest supersymmetric particle and be produced at LEP. 
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Neither e + e collider data [15| nor pp data from the Tevatron (T^j preclude this possi- 



bility. In contrast with Ref. [17 here we focus on an effective model where the breaking 
of R-parity is introduced through an explicit bilinear term in the superpotential. This is 
substantially simpler than the full majoron version of the model considered previously. 
In fact, this bilinear model is not only especially simple theoretically, also its phenomeno- 
logical implications in collider physics can be studied in a very systematic way. The 
bilinear model constitutes the simplest R-parity breaking model |27j consistent with ra- 
diative electroweak symmetry breaking, very much the same way as the minimal R-parity 



conserving supergravity models with universal soft SUSY breaking terms [IS|, MSUGRA, 
for short. As mentioned it also provides an attractive joint explanation of the present 
neutrino anomalies HI. 



In order to discuss stop decays we also refine the work presented in Ref. |L9|, |20], [21], ^2 
by giving, for the first time, an exact numerical calculation for the FCNC process t — > c\\. 
We also compare the results obtained this way with those one gets by adopting the usual 
one-step or leading logarithm approximation in the RGE's. In contrast with the R-parity 
conserving model such an approximation would be rather poor for our purposes, since 
we will be interested in comparing FCNC with R-parity violating stop decay modes (see 
section 5). Moreover, in contrast to ref. |I7| , where the magnitude of the stop - charm 
- neutralino coupling was a phenomenological parameter, here we assume a minimal 
supergravity scheme with universality of soft terms at the unification scale in which this 
coupling is induced radiatively. As we will see this has important phenomenological 
implications as for the behaviour of the R-parity violating stop decays with respect to 
tan/3. We calculate its magnitude using a set of RGE's in which the running of the 
Yukawa couplings and soft breaking terms is taken into account. Here we also provide 
the analysis of the relationship of the R-parity violating stop decays with the magnitude 
of the tau neutrino mass. Motivated by the simplest oscillation interpretation of the 
Super-Kamiokande atmospheric neutrino data, we also generalize the treatment of the 
R-parity violating decays by explicitly considering the case of light tau-neutrino masses, 
not previously discussed. 

For definiteness and simplicity we focus on supersymmetric models where the breaking 
of R-parity is parametrized explicitly through a bilinear superpotential term of the type 
£H U |[23|| . The stop can have new decay modes such as 

h^br (4) 

due to mixing between charged leptons and charginos. We show that this decay may be 
dominant or at least comparable to the ordinary R-parity conserving mode 

ti^cxl, (5) 

where Xi denotes the lightest neutralino. 

The paper is organized as follows. The model and an analytical analysis of the tree- 
level tau-neutrino in terms of SUGRA parameters is described in section |2|. The mass 
matrices are given in section |3] while in section 4 we present the top squark decay widths in 
the minimal supergravity model with universal soft SUSY breaking terms fl8| , MSUGRA, 



for short. The relevant Feynman rules and the squark decay widths and branching ratios 
are calculated in appendix |A|. They are studied numerically in section 5 and we present 
our conclusions in section 6. 



2 



2 The Model 



The supersymmetric Lagrangian is specified by the superpotential W given by 



w = e ab hfiosu/H* • h)\Q';i),ii';~ • „//,/.;' // ; 



fazrb 



[ b 



(6) 



where i,j = 1,2,3 are generation indices, a,b — 1,2 are £77(2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with e\ 2 = 1. The symbol "hat" over each letter 
indicates a superfield, with Qi, Li, Hi, and H 2 being SU{2) doublets with hypercharges 
—1, —1, and 1 respectively, and U, D, and R being SU{2) singlets with hypercharges 
— |, |, and 2 respectively. The couplings h v , hn and h E are 3x3 Yukawa matrices, and 
/1 and €i are parameters with units of mass. 

Supersymmetry breaking is parametrized by the standard set of soft supersymmetry 
breaking terms 

V so f t = MfQTQ* + MfU*Uj + MfD*Dj + M$ 2 L?L$ + M^ 2 R*Rj 

- ±M 3 A 3 A 3 + \M\ 2 \ 2 + iM'AjAi + h.c' 

+e ab {AZh%Q«U J H b 2 + A%h%Q h fi 3 Hl + A^p^F? 

-BilHIHI + B&LIH^ , (7) 

For definiteness and simplicity we assume only R-parity Violation (RPV) in the third 
generation, neglecting the effects of RPV on the two first families, adopting the superpo- 
tential @§@ 



W = h t Q 3 U 3 H 2 + h h Q 3 D 3 H x + h T L 3 R 3 H x - ^H X H 2 + e 3 L 3 H 2 (8) 

to describe the R-Parity violating violating decay mode t\ — > br. In this case we will 
omit the labels i, j in the soft breaking terms given above. Note that the bilinear term e 3 
can not be rotated away, since the rotation that eliminates it reintroduces an R-Parity 
violating trilinear term, as well as a sneutrino vacuum expectation value. Notice that, in 
contrast with ref. || where the doublet sneutrino vev in the bilinear model is much more 
loosely constrained, in this case it is not subject to constraints from astrophysics. 

Note, in contrast, that in order to describe Flavour Changing Neutral Current (FCNC) 
effects such as the R-Parity conserving process ti — > cx\ we need the three generations 
of quarks. 

The above model can be described in various equivalent bases, for example 

1. one in which bilinear term and sneutrino vev are non-zero, e 3 7^ and i> 3 7^ P7|, 0] 

2. one in which trilinear term[] and sneutrino vev are non-zero, A 3 7 7^ and v^ 1 7^ [2£ 

3. the vev-less basis in which e 3 7/ and X^ n are non-zero but t> 3 7/ = [^, [3(| 



the one generation case there is only one trilinear RPV term in the superpotential written in our 
notation as \3Q3D3L3 
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where the R-parity violating parameters can be expressed in terms of dimension-less basis- 
independent alignment parameters sin£, sin£' and sin£" [f2(], |3T|] (X = I, II or III) as 
follows: 

sme = 3 , ; = 4- = -V 9 

sin^ = ^ ^ = f = 4 (10) 

u d a b u d '% 

where 



K = yjhf + A? 2 // = V^ 2 + ^ = V^ 2 + «f *, X = /, J/, or III 

(12) 



Note that, in the notation of eqs. (p|)-(pT|), the parameters e 3 and /x appearing in 
eq. (§) should bear the superscript I. 

Of these parameters only two are independent because they satisfy 

sin £" = cos £' sin £ — sin £' cos £ (13) 

In the limit when the R-parity violating parameters vanish one recovers the MSSM. 
From now on we will work in the A3 = 0-basis, unless otherwise stated. As a result we will 
omit the label I in all the parameters associated with this basis. We also will drop out the 
prime in hb- One of the advantages in working in this basis is that the RGE's evolution 
does not induce the trilinear R-parity violating terms neither in the superpotential nor in 
the scalar potential if at the beginning we impose universality . 

It is convenient to introduce the following notation in spherical coordinates for the 
vacuum expectation values (vev): 

Vd = v sin 9 cos (5 
v u = v sin 9 sin j3 

t>3 = vcos9 (14) 

which preserves the standard MSSM definition tan/5 = v u /vd- In the MSSM limit, where 
£3 = V3 = 0, the angle 9 is equal to n/2. This makes sense in the A3 = 0-basis where the 
usual MSSM relation ^ 

fn b = —^hvd, (15) 

holds. 

In this model the presence of RPV induces a mass for the tau-neutrino at the tree 
level JO], [7], as well as radiative masses to the the v e and . As already mentioned it is 
sufficient for our present discussion of stop decays to keep only the tau-neutrino. 

In order to study the v T mass it is convenient to have an analytical expression for 
m Ur in this limit. The tree level tau-neutrino mass may be expressed as || [5], |26| - [32] 



m, 



(g 2M ' + W (16) 
4MM'fi' 2 -2(g 2 M' + g' 2 M) fx' v u v' d cos £ d ^ K J 
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in terms of basis-independent parameters fjf, v' d and sin£ defined in Eqs ([12]) and @. 
The second term in the denominator may be neglected if M, /i > mz, as often happens in 
minimal supergravity models with universal soft SUSY breaking terms WE\ . Thus one may 



obtain an estimate of the neutrino mass by keeping only the first term in the denominator. 



m. 



i2 



2M 



sin f 



(17) 



where we have used M' = Mg' 2 /g 2 . For sin£ « 1 one can easily check that m„ T could be 



as large as the experimental upper bound of 18 Mev [33]. However in MSUGRA models 
one may obtain naturally small sin £ values, calculable from the RGE evolution from the 
unification scale down to the weak scale. Indeed, using the minimization equations sin£ 
can be written in terms Am 2 



mjj — mi and AB 



B 3 - B [34] as 



sin£ 



/ Am 2 v 2f x'AB\ 
cos 4 sm 4 cos £ — 7) — I - 



m'lo 



v, m 



,2 



One may give a simplified approximate analytical expression for the tau-neutrino mass in 
this model by solving the renormalization group equations for the soft mass parameters 



m Hl , m Lz , B, and B3 in the one-step approximation. This gives [34| 



sin£ 



Am 2 



— cos £' sin £' cos £ h 2 

— cos £' sin £' cos £ 



^ + + M 2 D + A 2 D 



m'lo 



8tt 2 



In 



Mu 



■in" 



In 



(7 



and 



sin£|A# ~ cos£' sin £' tan 



1? 



m'?o 



In 



Mi 



u 



m t 



(19) 



(20) 



where we have denoted by the symbols sin£|Am 2 and sin£|As the two terms contributing 
to sin£ in eq. (IH). Using these expressions and assuming no strong cancellation between 
these terms one finds that the minimum neutrino mass is controlled by the sin£|Am 2 - As 
a result one finds, 



m, 



g 2 m\ 



(sin 2 £7i 



In 



Mr 



f21) 



M v ~ ~' \~6V m t 
The above approximate analytical form of the tau-neutrino mass is useful, as we will 



see later (e.g. eq. (|42|)) in order to display explicitly the degree of correlation between 
the R-parity violating decays, such as t\ — > br, with the tau-neutrino mass. 

The minimum value for sin£'/t& is determined by the value sin£' and that of tan/3. 
For sin £' ~ 1 and relatively small tan (3 so that h t is perturbative, one has 



77v Z lOKeV 



(22) 



for M ~ 1 TeV. In order to get smaller v T masses one needs to suppress sin 2 £' additionally, 
for example to reach one electron-volt the required R-parity violating parameters are 
given in Table [I|. These order-of-magnitude estimates are given in terms of the basis- 
independent angles £ and and in the relevant parameters for the three bases defined 
before. 
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basis-independent 


Basis I: A| = 


Basis II: e{/ = 


Basis III: 7 = 


sin£ 


sin f 


e 3 (GeV) 


vs(GeV) 




t^(GeV) 




ef'(GeV) 


1(T 5 


lCT 4 


icr 2 


IO" 4 


1 


icr 2 


1 


1(T 3 


IO" 4 


1(T 3 


IO" 4 


10~ 3 


IO" 2 



Table 1: Estimated magnitude of R-parity violating parameters required for a tau-neutrino 
mass in the eV range, without requiring cancellation in sin£ in the three bases defined before. 

Note that whenever the parameter has two values, the first correspond to tan/? = 2 
(the lower perturbativity limit) and the second to tan/3 = 35. In Table |l[ sin£ was 
estimated from eq. (|TT| ) and sin£' from eq. (pl|). 

Note that the RGE-induced suppression depends basically in the h% factor in eq. Q19D 
which is ~ (~ 1) for small (large) tan/5. As a result the bigger the value of tan/3, 
the smaller sin£' will have to be for a fixed tau neutrino mass. The RPV parameters in 
the several bases were estimated from Eqs. (PHTTI) and (O). 

In eq. ([H]) we have neglected AB contribution with respect to the one coming from 
Am 2 . It is possible, however, that the AB term may be sizeable. In the large AB case 
then it may cancel the Am 2 contribution in sin£, leading to an additionally suppressed 
neutrino mass. As we will see, however, in SUGRA models with universal soft terms at 
the unification scale (eSUGRA for short) we do not need any substantial cancellation in 
order to obtain v T masses below the electron-volt scale. 



3 Squark Mass Matrices 

The up and down-type squark mass matrices of our model have already been given previ- 
ously in Ref. |27j]. Here we generalize those to the three-generation case. The mass matrix 
of the up squark sector follows from the quadratic terms in the scalar potential 



V, 



quadratic 



fit T,t 



U 



U 



R 



Ml 



+ ■ 



(23) 



given by 



Ml 



1 

V2 



M 2 Q + \v 2 u huhu ] + A UL 



Ml + \vlhjhu + A UR 



(24) 



where A UL = \(g 2 - \g' 2 ){v 2 d - v 2 u + and A VR = y 2 (vj - v 2 u + vf)l are the splitting 
in the squark mass spectrum produced by electro-weak symmetry breaking, and A Vi - 



A l ^jh % jj. The eigenvalues of Ml are 



diag{m il ,m fi2 , . . . ,m fis } 



UL 



UR 



Ml 



1 UL 
1 UR 



(25) 



This way the six weak-eigenstate fields Un and um (i = 1, 2, 3) combine into six up-type 
mass eigenstate squarks Uk as follows: un = T^ k L u k = r^Wfc, u iR = T^ R u k = T*,f R Uk- 
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For completeness, we also give the mass matrix of the down squark sector. The 
quadratic scalar potential includes 



quadratic 



d* L ~d* R 



M 2 D 



d R 



+ ■ 



(26) 



given by 



M 2 D 



m q + \v 2 d h D h D ] + A DL -^v d A h D - -^fiv u h D 



\h t 



-^VdA^ - ^v u h D ] Ml + lv 2 d h D <<h D + A DR 



(27) 



where A DL = Adr = -±g' 2 (vj-vl+v!)l, and A Di - ee A%h%. 

The eigenvalues of Ml are 



diB l g{m 3i ,m 32 ,...,m d6 } 



r 



DL ^ DR 



Ml 



1 DL 

F f 

1 DR 



(28) 



One is left with six mass-eigenstate down squarks fields c4 related to dn and dm fields as 
follows: d iL = Tf L d k = T* D k id k} d iR = rf R d k = T* v k R d k . 

For the Higgs-slepton part of the quadratic scalar potential in the one generation case 
of the Bilinear R-parity Violating (BRpV) model, see refs. |35| and |23] . 

Of particular interest to us is the chargino/tau mass matrix. For our present purposes 
it is sufficient to have the form of this matrix for one generation, which is given by 



M, 



c 



M j= 2 gv 2 







(29) 



This form is common to all models with spontaneous breaking of R-parity, as well as in 
the simplest truncation of these models provided by the BRpV model considered here. 
We note that the chargino sector decouples from the tau sector in the limit €3 = v 3 = 0. 
As in the MSSM, the chargino mass matrix is diagonalized by two rotation matrices U 
and V 

m-± 
U*M C V- 1 = 1 m-± 
2 m T 



(30) 



The lightest eigenstate of this mass matrix must be the tau lepton (t^) and so the mass 
is constrained to be 1. 777051^26 GeV. To obtain this the tau Yukawa coupling becomes 



a function of the parameters in the mass matrix, and the full expression is given in [35]. 
The composition of the tau is given by 



<R 



where ip +T 



Hl 



and ip~ 



T, 



U^j (31) 



(-i\-,Hl 



r° ). The two-component Weyl 



spinors r R ~ and r° + are weak eigenstates, while t r and t~[ are the mass eigenstates. It 
follows easily from eq. ( |3"0"D that the matrix McMq is diagonalized by U and the matrix 
M^Mq is diagonalized by V. 
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The soft SUSY breaking parameters at the electroweak scale needed for the evalu- 
ation of the mass matrices and couplings are calculated by solving the renormalization 
group equations (RGE's) of the MSSM and imposing the radiative electroweak symmetry 
breaking condition. From the measured quark masses, CKM matrix elements and tan j3 
we first solve one-loop RGE's for the gauge and Yukawa couplings to calculate their cor- 
responding values at the unification scale. Assuming now universal soft supersymmetry 
breaking boundary conditions, we evolve downwards the RGE's for all MSSM parame- 
ters, including full three-generation mixing in the RGE's for Yukawa coupling constants, 
as well as soft SUSY breaking parameters. Next, we evaluate the Higgs potential at the 
m t scale including the one-loop corrections induced by the Yukawa coupling constants of 
the third generation. The radiative electroweak symmetry breaking requirement fixes the 
magnitude of the SUSY Higgs mass parameter \x and the soft SUSY breaking parameters 
B and B3. Notice that due to the third minimization condition one can solve for V3 as 
a function of €3. At this point, all RPV parameters at the electroweak scale are deter- 
mined as functions of the input parameters (tan/3, m , A , m^, sign(/i), e 3 ). Iteration is 
required because \i and e 3 are inputs to evaluate the loop-corrected minimum. Having 
determined all parameters at the electroweak scale, we obtain the masses and the mixings 
of all the SUSY particles by diagonalizing the corresponding mass matrices. At this stage 
we also choose 63 in order to get a sufficiently light tau-neutrino. 

We scan the soft SUSY breaking parameter space in the range 



100 GeV < 



m 


< 700 GeV 




< 400 GeV 


14)1 


< 1000 GeV, 


M 


< 200 GeV 


tan p 


< 60 



(32) 



the previous range on tan/3 guarantee that both ht and hb will be perturbative. For 
the CKM matrix, we use the Particle Data Group convention []3"E| ], taking V us = 0.2205, 
V c b = 0.041, IV^/V^I = 0.08 and neglecting CP violation, i.e. 5 = 0. Notice that here we 



scan over a much larger range for epsilon than used in ref. |L7 . 

The resulting region of lightest stop and chargino masses is displayed in Fig. [I]. Ne- 
glecting the three-body decays, we find that in Region I of the m^-m^r plane, BR(t\ — > 

cxl) + BR{h -> br) « 1. In Region II BR(t x -> br) + BR(t x -> bxf) « 1 (i=l,2). In 
Region III BR{t t -> brf + BRfa -> b + xt) + BR(t 1 -> tv T ) w 1 (i=l,2), while in region 
IV BR{t x -^br) + BR(t! -> bxf) + BR{t t -+tu T ) + BRfa -> t xj) « 1 (j = 1, • • • , 4). 
Note that in each region the exact equality to 1 is reached when the FCNC processes are 
fully included. 

In Appendix |A| we give the Feynman rules for all vertices involving squarks, quarks, 
charginos and neutralinos, as well as the two-body squark decay-widths, for squarks of all 
three generations. These equations reduce to the expressions found in Ref. |37| provided 
one identifies r^ 3 L = cos6^ and Tff R = sin%. They also generalize the results for the 
BRpV model to the three-generation case. 
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Figure 1: Kinematical regions in the mr -m~+ plane. From left to right: Region I mr < 

1 X l 1 

m-+ + m b ; Region II m~+ + m b < m?, < mt] Region III mt < mr < m^o + mt; and region IV 

> m^o + m t 

4 Lightest Stop Two-Body Decays in MSUGRA 



In an R-parity conserving supergravity theory the main t\ decay channel expected in 
region I of Fig. [I] is the loop-induced and flavour-changing t\ — > ex? [0, [2^, |2Tfl . 
As is well-known, the FCNC processes in the MSSM in general involve a very large 
number of input parameters. For this reason, following common practice, we prefer to 
perform the phenomenological study of flavour changing processes in the framework of a 
supergravity theory with universal supersymmetry breaking. The simplest description of 
FCNC processes in R-parity conserving minimal SUGRA models uses the so-called one- 
step approximation. Here we start by reproducing the standard calculation for t± — > c\\ 
as in [pl|. To do this consider only the effect of the third generation Yukawa coupling. 
From our general eq. ([A. 10|) we have for t\ = u\ 



r(ti 



with 



ex!) 



g 1 



57T 



(T 



UL13 



| sin BwN'u 



i 2 c ;n 2 a 



N' 



12 



COs9 W 



1 - 



m-o 

2 

ra- 
ti 



UL13 



A L cos 9r — A R sin 9~ 



2 

ra- 
th 



In the one-step approximation Ax,, Ar are given by 

tu 



(Mi 



lf)23 



{Mi 



123 



A 



(M? 



167T 2 

tu 



K cb K tb h 2 b (M 2 Q + M 2 D + ra 2 Hi + A, 
K cb K tb h 2 b ra t {A b + -A t ) 



with t 



u 



16vr 2 — ■ 2 - 

ln(Mc/m 4 ). So, in the one-step approximation we have 



r(ti 



cx?) 



Fh b (5 m 2 cos 0+ — 5a sin 0j 



(tan^iVn + 3iVi2) 



(33) 
(34) 

(35) 
(36) 

(37) 
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2 / \ 2 

where the pre-factor F = (^^K^Ktb) ~ 6 x 10~ 7 and the parameter S m 2 is given by 

5 m 2 = £ El i „ i (38) 

is basically independent of the initial conditions due to the mo dependence both in the numerator 
as in the denominator and 

_ MA + jA,) 
m~ — m~ 

Note however, that the one-step approximation includes only the third generation Yukawa 
couplings and neglects the running of the soft breaking terms ||, Such an approxi- 



mation is rather poor for our purposes, since we will be interested in comparing with R-parity 
violating decay modes (see the next section). In order to have an accurate calculation of the 
respective branching ratios we need to go beyond the one-step approximation. We therefore 
use a exact numerical calculation for the FCNC process t — ► c\i in which the running of the 
Yukawa couplings and soft breaking terms is taken into account. First we have checked that 
indeed the effect of the Yukawas from the two first generations is negligible. However the same 
is not true for the running of the soft breaking terms. As can be seen from Figure |2| the range 
of variation that we obtain from the numerical solution is 

r(ti -» ex?) ~ (10" 16 - l(T 6 )GeV (40) 

depending on the assumed value of m^j and tan/3. In this figure we have compared the decay 
width obtained from eq. ( |A.10| ) with the approximate formulae in eq. (|33| ) for two fixed values 
of mi/2, tan/? and taking Aq = 0. The approximate formulae only reproduce well the numerical 
result for the academic case of no SUSY breaking gaugino mass, m\/2 = 0. For the more realistic 
case mx/2 > 100 GeV, the exact solution is usually one decade smaller than the approximate one. 
In the one-step approximation Y(t\ — > cx±) can be arbitrarily small if the two terms 5 m 2 cos0£ 
and 5a sin 9^ in eq. (37) cancel. This behaviour can be illustrated in Figure |2| by the dashed line 



labelled 358, which corresponds to mo = 358 GeV. One sees clearly that while the approximate 
solution goes to zero, the numerical one reaches a minimum value around 10 -11 GeV. The wrong 
behaviour of the approximate solution indicates that the 5a depends strongly on the scale. For 
example, the RGE for At, is very sensitive on m±/2 an d tan/3 and in the one-step approximation 
there is no explicit dependence on m^, which is crucial. Both solutions increase with tan/3, as 
expected by the bottom Yukawa dependence explicit in eq. (|37|) and remain practically constant 
for large m^ values. 



5 Two-Body Decays of the Lightest Stop: the R- 
parity violating case 

In contrast to the case of an R-parity conserving supergravity theory, in our broken R-parity 
model one can have a competing R-parity violating stop decay mode in region I of Fig. |l|. From 
eq. flA.llD with r = F 3 + one can easily compute the R-parity violating stop decay width t\ — ► br, 

g 2 A 1//2 (m| ,m|,m 2 ) 
r(ti-*&T) = 1 3 {-4L^ 2 /i b ce t -(V3* 2 /i t S0 t - - V 3 * 1 ce t -)m f) m r 

h 

+m 2 h t s h - V 3 \c 6i ) 2 + U^hlc 2 e ,](ml -m 2 b - m 2 T )} (41) 
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Figure 2: Comparison between the exact numerical calculation (ordinate) and the one-step 
approximation (abscissa) for the t\ — ► ex? decay width for various values of tan/3 and m 1( / 2 
with = and tjiq varying in the indicated range. The dotted left diagonal line would signify 
the equality between the estimates, while the right diagonal line would indicate one order of 
magnitude difference. Results of both estimates indicated in the lower right legend. More 
details are found in the text. 



which coincides with the result found in Ref. [17]. In [38] it was shown that, except for U32 
which determines the SU(2)-conserving mixing of the Higgsino with the left-handed r, all other 
mixing matrix elements V%i and U^i are proportional to v^ 1 and therefore to the tau-neutrino 
mass. Neglecting these terms we have from eq. (UT 



r(ti 



br) 



m b ,m; 



167rm? 

ti 



noting that, to a good approximation, 



• 2 t'u2 2 / 2 
■sin ih b c ei {m k 



ml 



(42) 



\U : 



32 



sin£'| 



(43) 



where £3 corresponds to the bilinear mass parameter in basis I. The lesson here is that the R- 
Parity violating decay rate T(ti — > br) is proportional to €3 or, equivalently, to sin 2 instead of 
sin 2 £, and thus not necessarily small, since it is not directly controlled by the neutrino mass. In 
other words, there can be cancellations in the latter but not in the R-parity violating branching 
ratio. 

The meaning of the factor sin 2 ^'/i 2 may also be seen in basis II, where 63 = 0. In this case 
v^ 1 is proportional to the tau-neutrino mass so that, as already mentioned, in this basis all the 
elements U%i and V%i are small [ p6f[ . Neglecting these terms, T(t\ — ► br) may be written directly 
from the interaction term t^bjiT^, which is induced by the trilinear term in the = 0-basis 
given in eq. (jlC|) as 

(e 3 /fi')h b = h b sin^ (44) 



A 3 



which is the factor in eq. (f42|). Note, however, that in our numerical calculation to be described 
in the next section we have used for T(ti — > br) the full expression given in eq. ( A.10D of the 
appendix. 
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In the next section we will determine the conditions under which the R-parity violating 
decay width T(ti — > br) can be dominant over the R-parity conserving ones, T(t\ — > cxi) and 
Tih^bxt). 

5.1 Region I 

Using the one-step approximation for T(ii — ► ex?) one finds 

(m-o \ 2 
i-^f 1 ] (45) 

Using the eq. ( f42[) and neglecting charm, tau and bottom masses we get 

r(ti -> 6t) sin 2 £' I m| J v ; 

Therefore — » ex?) will start to compete with T(t\ —* br) from sin£' < 5 x 10~ 3 (10~ 4 ) for 
tan/3 large (small). In Fig. || we compare BR{ti — > cxi) [calculated numerically from their 
exact formula in ( |A.10 )1 with BR(t\ — > br) within the restricted region of the m^-m^o plane 



where only those two decay modes are open. We consider different m UT values (these correspond 
to relatively small values of the R-parity parameters | €3 1 , 1 1?3 1 < 1 GeV). We vary the MSSM 
parameters randomly obeying the condition < nn~± + nib an d depict the corresponding 
region in light grey. The upper-left triangular region is defined by kinematics and corresponds 
to mi x < m^o + m c , so that BR{t\ — ► br) = 100 %. The lower-right grey corresponds to 
> m^o + m c when the sampling is done over the region defined by eq. (j3^). One notices 
from Fig. || that in the central region the dominant stop decay mode is t% — > 6 r with branching 
ratio BR(t\ — > br) > 0.9. The dotted lines in the light grey region indicate maximum v T mass 
values obtained in the scan. In the calculation of the v T mass, we have allowed only up to one 
order of magnitude of cancellation between the two terms which contribute to sin£. Therefore 
if the lightest stop only decays into the two modes considered here, the processes t\ — > 6r, will 
be important even for the case of very light tau-neutrino masses. 

We note however that we can use the limits obtained from leptoquark searches []39| in order to 
derive limits on the top-squark for our R-parity violating case. In particular, if BR(t± — > br) = 1 
stop masses less than 99 GeV are excluded at 95% of CL., under the assumption that the three- 
body decays of the stops are negligible. Therefore, the dark region in Fig. |3] would be ruled out. 
In ref. [Ji(T| we have determined the corresponding restrictions on the SUGRA parameter space. 

The dependence on the tau neutrino mass may be seen in Fig. || where the role played by 
tan (3 is manifest. In this figure we have shown BR(t± — > br) as function of the lighter stop mass 
for tau-neutrino mass in the sub-eV range, indicated by the simplest oscillation interpretation 
of the Super-Kamiokande atmospheric neutrino data. We have obtained such tau-neutrino 
mass values numerically, allowing only one decade of cancellation between the two terms that 
contribute to sin£ in eq. (|l8|). The degree of suppression for e^/fi obtained numerically agrees 
very well with the expectations from the approximate formula for the minimal tau-neutrino 
mass m eq. (H). In contrast with Ref. @, 

in our case BR(t\ — > br) decreases with tan (5. The 



reason for this difference is that here we take into account the fact that the mixing parameter 



Fuli3 obtained from the RGE depends on h% in eq. (|33D , while in ref. [17] was simply regarded 
as a phenomenological input parameter (called 5 there). 

The message from this subsection is that in our SUGRA R-parity violating model the R- 
Parity violating decay mode t\ — > br can very easily dominate the R-Parity conserving decay 
mode t\ — > cxi, even for very small neutrino masses. 
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Figure 3: Regions where the t\ — ► br decay branching ratio exceeds 90% in the m^-m^o 
plane for different m Vr values. The MSSM parameters are randomly varied as indicated in the 
text under the restriction mx < m~± + m^. The upper-left triangular region corresponds to 

< m^o + m c so that only the t\ — ► b r decay channel is open. The lower-right unshaded 
region corresponds to > m-+ + m\,. 
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Figure 4: BR{t\ — > br) as function of the lighter stop mass for tau-neutrino mass in the sub-eV 
range and two different values of tan/3 and e^/fj,. This prediction is natural in the sense that 
we have allowed only up to one order of magnitude of cancellation between the two terms that 
contribute to sin£. 
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Figure 5: Contours of BR(t\ — > br) > BR(ti — > bxt) i n the mr -m-+ plane for |v 3 | < 10 GeV. 
Three different maximum values for | €3 1 are considered: \e^\ < 40 GeV (dot-dash), | e 3 1 < 60 
GeV (dots), and | €3 1 < 80 GeV (dashes). The region where mi < + m~+ corresponds to the 
previously studied Region I. 

5.2 Region II 

In region II the R-Parity conserving decay mode t\ — * bxf is open (but not t\ — ► tv), and 
competes with the R-Parity violating mode t\ — > br. Replacing the subindex 3 by 1 on the 
diagonalization matrices U and V in eq. (|4~l]) we get the corresponding expression for T(ti — ► 
bxt)- in order to get an approximate expression for the ratio of the two main decay rates in 
this region, we note that in MSUGRA the lightest chargino is usually gaugino-like, implying 
that ~ 1. In addition, the lightest stop is usually right-handed, hence sin 2 % > cos 2 %. This 
way we find 

T(r) _ r(t~i-6r) ^ sin 2 fK* cos 2 6~ t R 



F(xT ) r(ti -» 6 xi ) (V^* cos % - Ff^i sin %) 2 + C/i|/i 2 cos 2 ^ 

where if is a kinematical factor depending on the lightest stop and chargino masses, and here 
we have defined ht,b = ht,b/9- The presence of the bottom quark Yukawa coupling indicates 
that large values of tan (3 are necessary to have large R-Parity violating branching ratios in this 
region. In fact, we have checked numerically with the exact expressions that in Region II (RII) 
V(t)/T(xx) ^ 1 only for large tan/? as we will see in the next figures. 

In Fig. H we show the regions in the m x ± — in^ plane where BR(i\ — > br) dominates 

over BR{t\ — > bxt)- In the upper-left region the decay mode t\ — > &Xi^ is n °t allowed and 
corresponds to Region I. Below and to the right of this zone, and above and to the left of three 
rising lines, lies region RII where r(r)/r(xi") > 1. The three lines correspond to \e^\ < 80 GeV 
(dashed), | £3 1 < 60 GeV (dotted), and \e$\ < 40 GeV (dot-dashed), respectively. The proximity 
to the upper-left zone indicates that the RPV decay dominates only close to the threshold where 
there is a high kinematical suppression of the R-parity-conserving one, through the factor K. 
Unlike the case of region I this requires large values of the RPV parameters. Note, moreover, 
that if the stops have a small mixing (cos#£ ~ 0), then r(r)/r(xi") <C 1 in RII. 

A simpler expression for the ratio of decay rates in eq. (EW) is obtained if we take V\\ ~ 1 
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Figure 6: Regions for T(ti — > br)/T(ti — > b\\) as a function of the tau-neutrino mass with the 
universality condition B = B% at the unification scale imposed at the 0.1% level as indicated. Its 
effect is to alter the maximum attainable tau-neutrino mass. The dot-dashed line corresponds 
to the case where AB = at the weak scale. 



and assume no kinematical suppression in eq. (f47| ) through the factor K: 

T(t) 



sin 2 hi . (48) 



Note that the presence of the parameter sin£' = e^/ jJ indicates that the R-Parity violating 
decay mode is not strictly proportional to the neutrino mass, but proportional to the BRpV 
parameter ef>. 

However generically we expect some correlation with the v T mass, especially in the case where 
the boundary conditions in the RGE are universal and there are no strong cancellations between 
two terms that contribute to sin£ as shown in Fig. ^. In this figure we plot the ratio T(t)/F(x~i) 
in RII as a function of the tau-neutrino mass. Both decay rates have been calculated numerically 
from the exact formulas. In this figure we have imposed both m 2 Hi = and B = B3 within 
0.1% at the GUT scale. Cancellation between the Am 2 and AB terms in the neutrino mass 
formula of eq. (|l^) are accepted only within 1 decade. As a reference we have drawn the line 
corresponding to AB = and Am 2 = Am 2 nin (Am 2 is negative and its magnitude is bounded 
from below by Am^ jn ) at the weak scale, which gives an idea of the value of the neutrino mass 
when there is no cancellation between the AB and Am 2 terms. 

We have imposed an upper bound on m UT at the collider experimental limit of the tau- 
neutrino mass, and have chosen fixed values of e^/fi = 1, 0.1, and 0.01. The allowed region 
for e-s/fi = 1 is above the dashed line. In the case of e^/fi = 0.1 (0.01) the allowed region lies 
enclosed between the solid (dotted) lines. The effect of tan /3 is to increase the ratio r(r)/r(x+): 
the minimum value of the ratio is obtained for tan (3 ~ 2 and the maximum corresponds to 
tan (5 ~ 60. The extreme values of tan (3 are dictated by perturbativity. 

A number of statistically less significant points appear outside the drawn regions in Fig. ^ 
and are not depicted. They correspond to points with — mj, — Tn~± < 10 GeV which appear 
above the diagonal line, and points with cos% < 0.1 which appear below the horizontal line 
corresponding to the lowest values of tan 0. In the last case, our approximation in eq. (|48|) does 
not work any more. On the other hand, eq. (|^) predicts very well the behavior of T(t) /T(xt) 
if cos Of > 0.1. For example for 63///= 1, or equivalently sin^' = l/y/2, we expect from eq. (|48| ) 
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Figure 7: Regions for T(ti — » &r)/T(ii — ► as a function of the tau-neutrino mass for 

different levels of cancellation between the two terms that contribute to the neutrino mass. 
We impose the universality condition m 2 H = M 2 at the unification scale, but B and B% are 
not universal. We take \es/[j\ = 1 (inside the dashed lines), | €3 / | = 0.1 (solid lines), and 
|e 3 //u| = 0.01 (dotted lines). 

a maximum value of order 1 for large tan j3 {hi, ~ 1) and a minimum value of order 10~ 3 for 
small tan/3 (/i 2 ~ 10 -3 ), and this is confirmed by Fig. ||. High values of the R-parity violating 
branching ratio for large €3 values are highly restricted for large tan f3. This can be understood 
as follows. In the case of e^/fi = 1 and tan/3 = 60 acceptable neutrino masses are obtained 
only if sin£ ~ 1. On the other hand, in this regime we find from eq. (|l~8| ) that the AB term is 
large because of the high value of tan f3, and that the Am 2 term is large because m 2 H ^ becomes 
negative and Am 2 = m 2 Hi — M\ grows in magnitude. This way, acceptable neutrino masses 
are achieved only with cancellation within more than one decade. In any case, we think that 
Fig. H is very conservative considering that in MSSM-SUGRA with unification of top-bottom- 
tau Yukawa couplings, the large value of tan (3 implies that a cancellation of four decades among 
vev's is needed. 

The width of the band in Fig. |6| reflects the degree of correlation between the ratio T(t\ — > 
br)/T{t\ — > bxf) and the neutrino mass under the mentioned conditions. Note that one would 
have an indirect measurement of the neutrino mass if this ratio were determined independently. 
The band will open to the left if one allows a stronger cancellation between the terms in AB 
and Am 2 . On the other hand it will open to the right if the universality between B and -B3 
is relaxed. This is shown in Fig. [?] where we plot the ratio T(t)/T(xi) m RH as a function 
of the tau-neutrino mass, but without imposing universality between B3 and B. If we accept 
cancellation within one decade between the AB and Am 2 terms , then the allowed region is at 
the right and below the corresponding dashed tilted line. If a larger degree of cancellation is 
accepted, the left boundary of the allowed region moves to the left as indicated in the figure, 
enhancing the R-parity violating channel. In addition if we accept only a decade of cancellation 
between the two terms that contribute to the tau-neutrino mass, then our approximate formula 
which predicts the minimum tau-neutrino mass in eq. ( pT[ ) works very well. 

In summary, in this subsection we have shown that even in region II, where the R-Parity 
conserving decay mode t\ — > bxf is also open, the R-Parity violating decay mode t\ — > br 
can be comparable to t\ — > bxf f° r large tan/3 and £3, and relatively close to the chargino 
production threshold. In general, this implies a large neutrino mass unless a cancellation is 
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Figure 8: Universality condition B = B% at the unification scale as a function of vlo- As tan/3 
increases, the allowed values of Aq are more constrained. 



accepted between the two terms contributing to the tree level neutrino mass. In addition, the 
non-universality of the B and B$ terms at the GUT scale does not increase appreciably the 
allowed parameter space, except at large tan (3. The main consequence of this non-universality 
is to restrict the allowed values of Aq at large tan (3. In the next subsection we study the effects 
introduced by the non- universality of m 2 Hi and m^ 3 . 

5.3 Effects of non— Universality 

We now study the effect of possible non-universality of soft-breaking SUSY parameters on our 
previous results. In particular, the non-universality between rrijj and m 2 L at the GUT scale. 
The Minimal SUGRA model, while highly predictive, rests upon a number of simplifying as- 
sumptions which do not necessarily hold in specific models due to the possible evolution of the 
physical parameters in the range from Mpi anc k to Mqut- Specifically, there are several models 
in the literature with non-universal soft SUSY breaking mass parameters at high scales. A recent 



survey can be found in [41|, where several models such as based on string theory, M-theory, and 
anomaly mediated supersymmetry are analyzed. For this reason we find interesting to explore 
here the effects of non- universal soft terms. 

The SUGRA spectra are typically found for given values of m^, mo, Aq, tan/3 and Sgn(/i). 
In our case we have in addition sin£' (or equivalently, 63). The value of V3 is determined by the 
previous parameters through the minimization conditions. In addition, a relation between Aq 
and the ratio B3/B at the GUT scale (which indicates the degree of universality) emerges. This 
relation can be seen in Fig. || for 63/// = 1 and the values tan/3 = 3, 40, and 60, for m/jj = M 2 . 
The relation becomes more restrictive as tan/3 is increased, starting from —1000 < Aq < 1000 
GeV allowed for tan (3 = 3, to a single Aq value compatible with unification for tan (3 = 60. 

Another way to enhance the R-parity violating channel, enlarging the band towards the left 
in Fig. |(| is by relaxing the universality between m\ and M 2 at the GUT scale. In Fig. || we 
plot the ratio m 2 Hi /Mf j at the weak scale as a function of the same ratio at the unification scale 
Mqut for tan (3 = 3. The shaded region is allowed, implying a maximum value for the ratio 
m 2 Hi /Mf j at the weak scale for a given value of the ratio at the GUT scale. We see from Fig. || 
that a relaxation of universality of 0.5% or more is enough to make (m'jj i /M 2 ) wea f : = 1 possible, 
meaning that smaller neutrino masses are attainable without having to rely on a cancellation 



17 




Figure 9: Comparison between the ratio m 2 H /M 2 at the weak and the unification scales for 
tan (3 = 3. Universality at the unification scale, m 2 H /M 2 = 1, implies a maximum value for this 
ratio at the weak scale. 



between the Am 2 and AB terms or small values of T(ii — > 6r). However as we increase tan/3 
the maximum value of m 2 Hi /Mf j decreases, and thus, the required non-universality between m,H 1 
and Ml at unification scale grows drastically. In Fig. 10 we show the ratio m 2 Hl /M 2 at the 
weak scale as a function of tan/3. We appreciate clearly the growing of |Am 2 | m j n with tan/3. 
We remind the reader that this kind of non-universality in the soft terms is not uncommon in 
string models p2| , or GUT models based on SU(5) |4j| or SO(10) Q] for example. There are 
in fact some SO (10) models for non- universality of the GUT scale scalar masses which naturally 
favour light neutrino mass [44]. 

The effect of non-universality it is also explored in Fig. [IT] where it is shown the relation 
between the neutrino mass and the parameter sin£ for e^/fi = 1. Two different bands are 
shown: one for tan/3 = 3 and T(t) /F(xf) = 2 x 10~ 3 , and a second one for tan/3 = 46 and 
T(t)/T(xi) = 0.4 ±0.2. The required degree of universality at the GUT scale is indicated inside 
the bands. For example, in order to have neutrino masses of the order of eV for tan (3 = 3, m 2 Hi 
needs to be at least 0.2% larger than M 2 . Similarly, for tan/3 = 46 we need a m 2 H twice as 
large as M 2 at the GUT scale in order to have neutrino masses of 1 eV. We stress the fact that 
for Fig. 11 we have conservatively accepted cancellation at the level of one order-of-magnitude 
only. 

In summary, the lesson to learn here is that non-universal soft SUSY breaking terms at the 
GUT scale have the potential of making it easier to reconcile sizeable R-Parity violating effects 
in the stop sector with very small neutrino masses, without resorting to cancellations. 



6 Conclusions 

We have studied the decays of the lightest top squark in SUGRA models with and without 
R-parity. We have improved the calculation for the decay i\ — > c\° by numerically solving 
the renormalization group equations (RGE's) of the MSSM including full generation mixing 
in the RGE's for Yukawa couplings as well as soft SUSY breaking parameters. The decay- 
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Figure 10: (m 2 Hi /Mf j ) evaluated at the weak scale as a function of tan/3. This ratio is always 
less than one and decreases with tan j3. 
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Figure 11: Minimum value of the tau-neutrino mass as a function of sin£ for different values 
of m^j/M^ at the GUT scale and two values of tan/3. The ratio £3/// is fixed to the indicated 
value, leading to a nearly constant value for T(i\ — > &T)/T(ii — > bxt)- Here we assume that 
the two terms contributing to the tau-neutrino mass cancel to within an order of magnitude. 



19 



width is in general one order of magnitude smaller than the one obtained in the usual one-step 
approximation. This result will therefore enlarge the regions of parameter space where the 
four-body decays of lightest stop dominate over the decay into a charm quark and the lightest 
neutralino. As a result it will affect the present experimental lower bound on the t\ mass even 
in the R-parity conserving case [p2j] . If R-parity breaks of course new decay modes appear 
and, as we have shown, they can be sizeable. In fact we have shown that the lightest stop 
can be the LSP, decaying with 100% rate into a bottom quark and a tau lepton. We have 
shown that the decay mode i\ —* br dominates over t\ — ► cx° even for neutrino masses in the 
range suggested by the simplest oscillation interpretation of the Super-Kamiokande atmospheric 
neutrino data. This result would have a strong impact on the top squark search strategies at 
LEP [45] and TEVATRON 46], where it is usually assumed that the t\ — > cx° decay mode 
is the main channel. In addition to the signal of two jets and two taus present when the two 
produced stops decays through the R-parity violating channel, one expects a plethora of exotic 
high-multiplicity fermion events arising from neutralino decay, since such decay can happen 
inside the detector even for the small neutrino masses in the range suggested by the Va to v T 
oscillation interpretation of the atmospheric neutrino anomaly H . We have also compared the 
decay t\ —* br with the R-parity and flavor conserving mode t\ — > b\ + and shown that the 
rate of the former can be comparable or even bigger than the latter if the tau-neutrino mass 
and tan/3 are large. However one may have a sizeable branching of t\ — ► br in the case of 
suppressed tree level neutrino mass as a result of strong cancellations between the two terms 
that contribute to sin£, or in some regions of parameter space of non- universal SUGRA models 
with (m'jj i /m'j j3 )GUT ^ 1. A detailed analysis of the detectability prospects of such related 
signatures at present and future accelerators lies outside of the scope of the present paper and 
it will be taken up elsewhere. 
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A Feynman Rules 



In this Appendix we derive the Feynman rules F^qiq^ (involving a neutralino/tau-neutrino, a 
quark, and a squark) and F^qiq' k (involving a chargino/tau, a quark, and a squark of different 
electric charge) in the case of three generations and RPV in the third generation. This is a 
generalization of the Feynman rules contained in [^?J , which are done for the R-Parity conserving 
MSSM and for one generation of quarks and squarks. 

Following [fig] we work in a quark interaction basis where <1l,r 
ur = u R (we denote q and q° the mass and current eigenstates respectively), as opposed to 
Ref. |4S] where a more general basis is used. In addition, we implement the notation <jl,r = Ql r 
for the interaction basis. 

The starting point is the following piece of the Lagrangian 



d L,R> u l 



Ku° L , and 



£-uuF° 



~9 u i 



V2 



sin OweuN'jx + 



cos 9w 



eu sin 2 0w)N'j2 
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0(7 



\^2m\Y sin (3 sin 9 



N'mUjAPrF*} 



sin 8\yeuN' ji + 



1 



cos 9w 



-ejj sin 2 6yy)N'* 



J2 



UiR 



in 



\/2mw sin (3 sin 8 



N'* M u jL }P L F° + h.c. 



(A.l) 



written in the quark interaction basis. The 5x5 matrix N' diagonalizes the neutralino/neutrino 
mass matrix in the (7, Z, H®, H®, v r) basis as defined in f35fl , with the index J = 1...5. The 3x3 
up-type quark mass matrix m ou is not diagonal, with the indexes i,j = 1, 2, 3. 

In order to write the above Lagrangian with mass eigenstates we use the basic relations 
mentioned before, in particular, v9 L = {K'Y^UjL, which implies that u® L = ujlK^ 1 . We need 
the following relations: 



U° iL UiL 
UiRUiR 



u iL (T* UR m u ) ki u k 
u iR T$ R u k 

u iR (T UL K*m u ) ki u k 



(A.2) 



where i, j = 1,2,3 label the quark flavours, k = 1...6 labels the squarks, and m u = diag {m u , m c , mt} 
is the diagonal up-type quark mass matrix. In this way, the Lagrangian in eq. (|A.l[) can be 
written as 



-gu^G^l + H^ R )P R - (V2G^ R - H$i)P L ]F»u k + h. 
where the different couplings are 



(A.3) 



G 



,jki 
OUL 



smOweuN'^ + —^~{\ ~ eu sin 2 9 W )N'*J (v UL K^ 



n jki 
^OUR 

zrjki 
n OUL 

jrjki 
n 0UR 



COS 9\y 2 

sin OweuN'ji + — (— eu sin 2 9w)^'j2 



cos 6w 



UR 



(A.4) 



N' jA {T UL tfhi 
N'* jA (r UR hu) ki 



and hu =diag (m u , m c , mt)/{y2mw sin /? sin (9). Graphically, the FjUiU k Feynman rules are 
given by 




-igmc&L + H Iur)Pl ~ (V2G$i R - H^ L )P R ] 



Kv^JSg + H*j»)P R - - Kij k l)P L ] 



-ig 



The analogous Feynman rules in the MSSM are obtained by replacing Ff —> by interpretin. 
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the matrix N as the usual 4x4 neutralino mixing matrix, and by setting 9 = ir/2 in the formula 
for the Yukawa couplings. 

Similarly, replacing all u(u) by d(d) in eq. (|A.1| ) and starting from 



£qq'F+ 



gd 

+gdi 



ni 



ou\ 



\^2m\Y sin j3 sin 9 
m ?j 

\^2m\y cos (3 sin 9 



Vj2UjR - VjiU iL 



PrFj 



U} 2 u jL P l F c j 



m 



\/2mw cos j3 sin 9 



Uj2djR — Ujidii 



+gu1 



771 



OU 



\^2m\y cos (3 sin 9 



V] 2 b jL P L Fj + h.c 



(A.5) 



we can obtain the complete Feynman rules for the neutralino/tau-neutrino and chargino/tau 
with quarks and squarks. The results, that complements the obtained in [48], are 
Neutralino-(d) quark-(d) squark 




ig[(V2G^ L + H^ R )P L - (V2G^ R - H^ L )P R ] 



T jki 



jki 
ODR 



T jki 



ODR 



The mixing matrices God and Hqd are defined as 



jki 



G 



ODL 

,jki 
ODR 



rrjki 

ODL 
Tj*jki 
n ODR 



sm9\yeDN'* 1 + 

sin 9wcdN' ji + 

N' j3 (T DL h D ) ki 
N'UT DR h D ) H 



1 



cos 9w 



(T 3D - e D sin 2 9 W )N'* 2 



1 DL 



1 



cos 9w 



(-eD sin 2 9w)N' j2 



DR 



(A.6) 



Chargino/tau-(d)quark-(u)squark 
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ig(-C-i)[(Gi k i - H^ R )P L - HtflP R ] 



rjki - 




-ig[(G*J ki - H* v f)P R - H*J ki P L ]C 



where C is the charge conjugation matrix (in spinor space) and the mixing matrices Gjj and 
Hu are defined as 

G$L = V A T ul, Hij k l^U* 2 {T UL h D ) k \ 

H$Tr = V; 2 {T UR huK) k \ (A.7) 
Chargino/tau-(u)quark-(d)squark 




ig(-C-i)[(G%i - H$ R )P L - HtflPn] 



jki 



*jki 



where the mixing matrices Gd and Hp are defined as 

U* x {T DL K^) k \ H** = V* 2 (T DL K^hu) k \ 
, 1)H = U* 2 (T DR h D K^ 



r jki 
DL 



In order to derive the decays widths we write, for example eq. ( A. 3 ) as 



(Ai 



Arfii* = 9u1(fu jH PR + hTP L )F% k + h.c 



iki - 



(A.9) 



The result is 



F(q k -» qi + Ff) 



g 2 X 1 / 2 (m 2 g k ,m 2 q .,m 2 F o) 
167rm? 

Qk 



+ (ihff) 2 + (fi ki ) 2 )(m, 



m qi m F o 



"q k m Qi ~ m F° 



(A.10) 
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16"7rm~ 

Ik 

{iff + (<l) 2 



where Q = U,D refers to q and 



f Q 

. jki 

l Q 
, jki 



m 



Ik 



m„ 



m 2 F ± 
j 



OQRJ 



VZU 0QR 



OQL 



QR 



n jki 



(A.ll) 

(A.12) 
(A.13) 
(A.14) 



with the G and H couplings defined earlier in this appendix. 
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